Mapping interictal discharges using intracranial EEG-fMRI to predict postsurgical outcomes.

Various subjective and objective methods have been proposed to identify which interictal epileptiform discharge (IED)-related EEG-fMRI results are more likely to delineate seizure generating tissue in patients with drug-resistant focal epilepsy for the purposes of surgical planning. In this intracranial EEG-fMRI study, we evaluated the utility of these methods to localize clinically relevant regions pre-operatively and compared the extent of resection of these areas to post-operative outcome. Seventy patients admitted for intracranial video-EEG monitoring were recruited for a simultaneous intracranial EEG-fMRI study. For all analyses of blood oxygen level-dependent responses associated with IEDs, an experienced epileptologist identified the most Clinically Relevant brain activation cluster using available clinical information. The Maximum cluster (the cluster with the highest z-score) was also identified for all analyses and assigned to one of three confidence levels (low, medium, or high) based on the difference of the peak z-scores between the Maximum and Second Maximum cluster (the cluster with the second highest peak z-value). The distance was measured and compared between the peak voxel of the aforementioned clusters and the electrode contacts where the interictal discharge and seizure onset were recorded. In patients who subsequently underwent epilepsy surgery, the spatial concordance between the aforementioned clusters and the area of resection was determined and compared to post-operative outcome. We evaluated 106 different IEDs in 70 patients. Both subjective (identification of the Clinically Relevant cluster) and objective (Maximum cluster much more significant than the second maximum cluster) methods of culling non-localizing EEG-fMRI activation maps increased the spatial concordance between these clusters and the corresponding IED or seizure onset zone contacts. However, only the objective methods of identifying medium and high confidence maps resulted in a significant association between resection of the peak voxel of the Maximum cluster and post-operative outcome. Resection of this area was associated with good post-operative outcomes but was not sufficient for seizure freedom. On the other hand, we found that failure to resect the medium and high confidence Maximum clusters was associated with a poor post-surgical outcome (negative predictive value = 1.0, sensitivity = 1.0). Objective methods to identify higher confidence EEG-fMRI results are needed to localize areas necessary for good post-operative outcomes. However, resection of the peak voxel within higher confidence Maximum clusters is not sufficient for good outcomes. Conversely, failure to resect the peak voxel in these clusters is associated with a poor post-surgical outcome.

In situ theranostic platform combining highly localized magnetic fluid hyperthermia, magnetic particle imaging, and thermometry in 3D.

Theranostic platforms, combining diagnostic and therapeutic approaches within one system, have garnered interest in augmenting invasive surgical, chemical, and ionizing interventions. Magnetic particle imaging (MPI) offers a quite recent alternative to established radiation-based diagnostic modalities with its versatile tracer material (superparamagnetic iron oxide nanoparticles, SPION). It also offers a bimodal theranostic framework that can combine tomographic imaging with therapeutic techniques using the very same SPION. Methods: We show the interleaved combination of MPI-based imaging, therapy (highly localized magnetic fluid hyperthermia (MFH)) and therapy safety control (MPI-based thermometry) within one theranostic platform in all three spatial dimensions using a commercial MPI system and a custom-made heating insert. The heating characteristics as well as theranostic applications of the platform were demonstrated by various phantom experiments using commercial SPION. Results: We have shown the feasibility of an MPI-MFH-based theranostic platform by demonstrating high spatial control of the therapeutic target, adequate MPI-based thermometry, and successful in situ interleaved MPI-MFH application. Conclusions: MPI-MFH-based theranostic platforms serve as valuable tools that enable the synergistic integration of diagnostic and therapeutic approaches. The transition into in vivo studies will be essential to further validate their potential, and it holds promising prospects for future advancements.

Increased interictal synchronicity of respiratory related brain pulsations in epilepsy.

Respiratory brain pulsations have recently been shown to drive electrophysiological brain activity in patients with epilepsy. Furthermore, functional neuroimaging indicates that respiratory brain pulsations have increased variability and amplitude in patients with epilepsy compared to healthy individuals. To determine whether the respiratory drive is altered in epilepsy, we compared respiratory brain pulsation synchronicity between healthy controls and patients. Whole brain fast functional magnetic resonance imaging was performed on 40 medicated patients with focal epilepsy, 20 drug-naïve patients and 102 healthy controls. Cerebrospinal fluid associated respiratory pulsations were used to generate individual whole brain respiratory synchronization maps, which were compared between groups. Finally, we analyzed the seizure frequency effect and diagnostic accuracy of the respiratory synchronization defect in epilepsy. Respiratory brain pulsations related to the verified fourth ventricle pulsations were significantly more synchronous in patients in frontal, periventricular and mid-temporal regions, while the seizure frequency correlated positively with synchronicity. The respiratory brain synchronicity had a good diagnostic accuracy (ROCAUC = 0.75) in discriminating controls from medicated patients. The elevated respiratory brain synchronicity in focal epilepsy suggests altered physiological effect of cerebrospinal fluid pulsations possibly linked to regional brain water dynamics involved with interictal brain physiology.

Aphasia recovery by language training using a brain-computer interface: a proof-of-concept study.

Aphasia, the impairment to understand or produce language, is a frequent disorder after stroke with devastating effects. Conventional speech and language therapy include each formal intervention for improving language and communication abilities. In the chronic stage after stroke, it is effective compared with no treatment, but its effect size is small. We present a new language training approach for the rehabilitation of patients with aphasia based on a brain-computer interface system. The approach exploits its capacity to provide feedback time-locked to a brain state. Thus, it implements the idea that reinforcing an appropriate language processing strategy may induce beneficial brain plasticity. In our approach, patients perform a simple auditory target word detection task whilst their EEG was recorded. The constant decoding of these signals by machine learning models generates an individual and immediate brain-state-dependent feedback. It indicates to patients how well they accomplish the task during a training session, even if they are unable to speak. Results obtained from a proof-of-concept study with 10 stroke patients with mild to severe chronic aphasia (age range: 38-76 years) are remarkable. First, we found that the high-intensity training (30 h, 4 days per week) was feasible, despite a high-word presentation speed and unfavourable stroke-induced EEG signal characteristics. Second, the training induced a sustained recovery of aphasia, which generalized to multiple language aspects beyond the trained task. Specifically, all tested language assessments (Aachen Aphasia Test, Snodgrass & Vanderwart, Communicative Activity Log) showed significant medium to large improvements between pre- and post-training, with a standardized mean difference of 0.63 obtained for the Aachen Aphasia Test, and five patients categorized as non-aphasic at post-training assessment. Third, our data show that these language improvements were accompanied neither by significant changes in attention skills nor non-linguistic skills. Investigating possible modes of action of this brain-computer interface-based language training, neuroimaging data (EEG and resting-state functional MRI) indicates a training-induced faster word processing, a strengthened language network and a rebalancing between the language- and default mode networks.

Holo-Hilbert spectral-based noise removal method for EEG high-frequency bands.

Simultaneous EEG-fMRI is a growing and promising field, as it has great potential to further our understanding of the spatiotemporal dynamics of brain function in health and disease. In particular, there is much interest in understanding the fMRI correlates of brain activity in the gamma band (> 30 Hz), as these frequencies are thought to be associated with cognitive processes involving perception, attention, and memory, as well as with disorders such as schizophrenia and autism. However, progress in this area has been limited due to issues such as MR-induced artifacts in EEG recordings, which seem to be more problematic for gamma frequencies. This paper presents a noise removal method for the gamma band of EEG that is based on the Holo-Hilbert spectral analysis (HHSA), but with a new implementation strategy. HHSA uses a nested empirical mode decomposition (EMD) to identify amplitude and frequency modulations (AM and FM, respectively) by averaging over frequencies with high and significant powers. Our method examines gamma band by applying two layers of EMD to the FM and AM components, removing components with very low power based on the power-instantaneous frequency spectrum, and subsequently reconstructs the denoised gamma-band signal from the remaining components. Simulations demonstrate that our proposed method efficiently reduces artifacts while preserving the original gamma signal which is especially critical for simultaneous EEG/fMRI studies.

The Evolution and Future of Spinal Drains for Thoracic Aortic Aneurysm Repair: A Review.

For decades, spinal drains for cerebrospinal fluid (CSF) pressure monitoring and drainage have been used as adjuncts to protect against spinal cord injury resulting from thoracic aortic aneurysm repair. There are many different approaches to placement and management of CSF drains, with no true consensus on best practice. Furthermore, the incidence of complications resulting from spinal drains largely has been stagnant. This review describes the history and rationale behind placement of CSF drains, explore various considerations, techniques, and equipment, and discuss potential considerations for developing more comprehensive protocols.

Topography-Related EEG-fMRI in Surgically Confirmed Epileptic Foci: A Comparison to Spike-Related EEG-fMRI in Clinical Practice.

EEG-fMRI has gained increasing importance in epilepsy pre-surgical diagnosis. However, 40-70% of EEG-fMRI recordings in patients lack interictal epileptiform discharges (IEDs) during the scan, which could be overcome by detecting matching topography maps. We tried to validate this method in clinical settings taking various electroclinical factors into consideration. Eleven patients who had undergone EEG-fMRI during pre-surgical evaluation for drug-resistant epilepsy and who had had clinical long-term video-EEG were studied. Spike-related blood oxygen level-dependent (BOLD) maps were created using IEDs occurring during the EEG-fMRI scan. Separate maps were then generated from IEDs marked on the clinical long-term EEG recordings, which were averaged to produce topographical IED maps and correlated with the EEGs recorded inside the scanner yielding a correlation coefficient time course. Epileptogenic zones were defined by an expert panel during pre-surgical evaluation and validated by an epilepsy surgery resulting in a good outcome. Both techniques' performance was evaluated according to factors including arousal during IED recording, IED topography and lateralization, lesion type, and localization. Topography-related EEG-fMRI yielded more specific results compared to the spike-related method. Superficial lesion location and ipsilateral IED seem to result in a higher concordance of BOLD maps. The polarity of BOLD responses may be lesion-dependent, and both positive and negative BOLD changes may be associated with the irritative zone. Topography-related EEG-fMRI may show improved specificity especially for superficial lesions producing ipsilateral spikes. This method can be used as an alternative either in the absence of spikes during the simultaneous EEG-fMRI acquisition or to sharpen a diffusely activated BOLD-map.

Trading off spatio-temporal properties in 3D high-speed fMRI using interleaved stack-of-spirals trajectories.

PURPOSE: Highly undersampled acquisitions have been proposed to push the limits of temporal resolution in functional MRI. This contribution is aimed at identifying parameter sets that let the user trade-off between ultra-high temporal resolution and spatial signal quality by varying the sampling densities. The proposed method maintains the synergies of a temporal resolution that enables direct filtering of physiological artifacts for highest statistical power, and 3D read-outs with optimal use of encoding capabilities of multi-coil arrays for efficient sampling and high signal-to-noise ratio (SNR). METHODS: One- to four-shot interleaved spherical stack-of-spiral trajectories with repetition times from 96 to 352 ms at a nominal resolution of 3 mm using different sampling densities were compared for image quality and temporal SNR (tSNR). The one- and three-shot trajectories were employed in a resting state study for functional characterization. RESULTS: Compared to a previously described single-shot trajectory, denser sampled trajectories of the same type are shown to be less prone to blurring and off-resonance vulnerability that appear in addition to the variable density artifacts of the point spread function. While the multi-shot trajectories lead to a decrease in tSNR efficiency, the high SNR due to the 3D read-out, combined with notable increases in image quality, leads to superior overall results of the three-shot interleaved stack of spirals. A resting state analysis of 15 subjects shows significantly improved functional sensitivity in areas of high off-resonance gradients. CONCLUSION: Mild variable-density sampling leads to excellent tSNR behavior and no increased off-resonance vulnerability, and is suggested unless maximum temporal resolution is sought.

Improving the sensitivity of spin-echo fMRI at 3T by highly accelerated acquisitions.

PURPOSE: Spin-echo (SE) functional MRI (fMRI) can be highly advantageous compared to gradient-echo (GE) fMRI with respect to magnetic field-inhomogeneity artifacts. However, at 3T, the majority of blood oxygenation level-dependent (BOLD) fMRI experiments are performed using T2∗ -weighted GE sequences because of their superior sensitivity compared to SE-fMRI. The presented SE implementation of a highly accelerated GE pulse sequence therefore aims to improve the sensitivity of SE-fMRI while profiting from a reduction of susceptibility-induced signal dropout. METHODS: Spin-echo MR encephalography (SE-MREG) is compared with the more conventionally used spin-echo echo-planar imaging (SE-EPI) and spin-echo simultaneous multislice (SE-SMS) at 3T in terms of capability to detect neuronal activations and resting-state functional connectivity. For activation analysis, healthy subjects underwent consecutive SE-MREG (pulse repetition time [TR] = 0.25 seconds), SE-SMS (TR = 1.3 seconds), and SE-EPI (TR = 4.4 seconds) scans in pseudorandomized order applied to a visual block design paradigm for generation of t-statistics maps. For the investigation of functional connectivity, additional resting-state data were acquired for 5 minutes and a seed-based correlation analysis using Stanford's FIND (Functional Imaging in Neuropsychiatric Disorders) atlas was performed. RESULTS: The increased sampling rate of SE-MREG relative to SE-SMS and SE-EPI improves the sensitivity to detect BOLD activation by 33% and 54%, respectively, and increases the capability to extract resting-state networks. Compared with a brain region that is not affected by magnetic field inhomogeneities, SE-MREG shows 2.5 times higher relative signal strength than GE-MREG in mesial temporal structures. CONCLUSION: SE-MREG offers a viable possibility for whole-brain fMRI with consideration of brain regions that are affected by strong susceptibility-induced magnetic field gradients.

The neuronal associations of respiratory-volume variability in the resting state.

The desire to enhance the sensitivity and specificity of resting-state (rs-fMRI) measures has prompted substantial recent research into removing noise components. Chief among contributions to noise in rs-fMRI are physiological processes, and the neuronal implications of respiratory-volume variability (RVT), a main rs-fMRI-relevant physiological process, is incompletely understood. The potential implications of RVT in modulating and being modulated by autonomic nervous regulation, has yet to be fully understood by the rs-fMRI community. In this work, we use high-density electroencephalography (EEG) along with simultaneously acquired RVT recordings to help address this question. We hypothesize that (1) there is a significant relationship between EEG and RVT in multiple EEG bands, and (2) that this relationship varies by brain region. Our results confirm our first hypothesis, although all brain regions are shown to be equally implicated in RVT-related EEG-signal fluctuations. The lag between RVT and EEG is consistent with previously reported values. However, an interesting finding is related to the polarity of the correlation between RVT and EEG. Our results reveal potentially two main regimes of EEG-RVT association, one in which EEG leads RVT with a positive association between the two, and one in which RVT leads EEG but with a negative association between the two. We propose that these two patterns can be interpreted differently in terms of the involvement of higher cognition. These results further suggest that treating RVT simply as noise is likely a questionable practice, and that more work is needed to avoid discarding cognitively relevant information when performing physiological correction rs-fMRI.

15 Years MR-encephalography.

OBJECTIVE: This review article gives an account of the development of the MR-encephalography (MREG) method, which started as a mere 'Gedankenexperiment' in 2005 and gradually developed into a method for ultrafast measurement of physiological activities in the brain. After going through different approaches covering k-space with radial, rosette, and concentric shell trajectories we have settled on a stack-of-spiral trajectory, which allows full brain coverage with (nominal) 3 mm isotropic resolution in 100 ms. The very high acceleration factor is facilitated by the near-isotropic k-space coverage, which allows high acceleration in all three spatial dimensions. METHODS: The methodological section covers the basic sequence design as well as recent advances in image reconstruction including the targeted reconstruction, which allows real-time feedback applications, and-most recently-the time-domain principal component reconstruction (tPCR), which applies a principal component analysis of the acquired time domain data as a sparsifying transformation to improve reconstruction speed as well as quality. APPLICATIONS: Although the BOLD-response is rather slow, the high speed acquisition of MREG allows separation of BOLD-effects from cardiac and breathing related pulsatility. The increased sensitivity enables direct detection of the dynamic variability of resting state networks as well as localization of single interictal events in epilepsy patients. A separate and highly intriguing application is aimed at the investigation of the glymphatic system by assessment of the spatiotemporal patterns of cardiac and breathing related pulsatility. DISCUSSION: MREG has been developed to push the speed limits of fMRI. Compared to multiband-EPI this allows considerably faster acquisition at the cost of reduced image quality and spatial resolution.

The variability of functional MRI brain signal increases in Alzheimer's disease at cardiorespiratory frequencies.

Biomarkers sensitive to prodromal or early pathophysiological changes in Alzheimer's disease (AD) symptoms could improve disease detection and enable timely interventions. Changes in brain hemodynamics may be associated with the main clinical AD symptoms. To test this possibility, we measured the variability of blood oxygen level-dependent (BOLD) signal in individuals from three independent datasets (totaling 80 AD patients and 90 controls). We detected a replicable increase in brain BOLD signal variability in the AD populations, which constituted a robust biomarker for clearly differentiating AD cases from controls. Fast BOLD scans showed that the elevated BOLD signal variability in AD arises mainly from cardiovascular brain pulsations. Manifesting in abnormal cerebral perfusion and cerebrospinal fluid convection, present observation presents a mechanism explaining earlier observations of impaired glymphatic clearance associated with AD in humans.

Corrigendum: Histological Correlates of Diffusion-Weighted Magnetic Resonance Microscopy in a Mouse Model of Mesial Temporal Lobe Epilepsy.

[This corrects the article DOI: 10.3389/fnins.2020.00543.].

Respiratory-related brain pulsations are increased in epilepsy-a two-centre functional MRI study.

Resting-state functional MRI has shown potential for detecting changes in cerebral blood oxygen level-dependent signal in patients with epilepsy, even in the absence of epileptiform activity. Furthermore, it has been suggested that coefficient of variation mapping of fast functional MRI signal may provide a powerful tool for the identification of intrinsic brain pulsations in neurological diseases such as dementia, stroke and epilepsy. In this study, we used fast functional MRI sequence (magnetic resonance encephalography) to acquire ten whole-brain images per second. We used the functional MRI data to compare physiological brain pulsations between healthy controls (n = 102) and patients with epilepsy (n = 33) and furthermore to drug-naive seizure patients (n = 9). Analyses were performed by calculating coefficient of variation and spectral power in full band and filtered sub-bands. Brain pulsations in the respiratory-related frequency sub-band (0.11-0.51 Hz) were significantly (P < 0.05) increased in patients with epilepsy, with an increase in both signal variance and power. At the individual level, over 80% of medicated and drug-naive seizure patients exhibited areas of abnormal brain signal power that correlated well with the known clinical diagnosis, while none of the controls showed signs of abnormality with the same threshold. The differences were most apparent in the basal brain structures, respiratory centres of brain stem, midbrain and temporal lobes. Notably, full-band, very low frequency (0.01-0.1 Hz) and cardiovascular (0.8-1.76 Hz) brain pulses showed no differences between groups. This study extends and confirms our previous results of abnormal fast functional MRI signal variance in epilepsy patients. Only respiratory-related brain pulsations were clearly increased with no changes in either physiological cardiorespiratory rates or head motion between the subjects. The regional alterations in brain pulsations suggest that mechanisms driving the cerebrospinal fluid homeostasis may be altered in epilepsy. Magnetic resonance encephalography has both increased sensitivity and high specificity for detecting the increased brain pulsations, particularly in times when other tools for locating epileptogenic areas remain inconclusive.

Analysis of accelerated 4D flow MRI in the murine aorta by radial acquisition and compressed sensing reconstruction.

Preclinical 4D flow MRI remains challenging and is restricted for parallel imaging acceleration due to the limited number of available receive channels. A radial acquisition with combined parallel imaging and temporal compressed sensing reconstruction was implemented to achieve accelerated preclinical 4D flow MRI. In order to increase the accuracy of the measured velocities, a quantitative evaluation of different temporal regularization weights for the compressed sensing reconstruction based on velocity instead of magnitude data is performed. A 3D radial retrospectively triggered phase contrast sequence with a combined parallel imaging and compressed sensing reconstruction with temporal regularization was developed. It was validated in a phantom and in vivo (C57BL/6 J mice), against an established fully sampled Cartesian sequence. Different undersampling factors (USFs [12, 15, 20, 30, 60]) were evaluated, and the effect of undersampling was analyzed in detail for magnitude and velocity data. Temporal regularization weights λ were evaluated for different USFs. Acceleration factors of up to 20 compared with full Nyquist sampling were achieved. The peak flow differences compared with the Cartesian measurement were the following: USF 12, 3.38%; USF 15, 4.68%; USF 20, 0.95%. The combination of 3D radial center-out trajectories and compressed sensing reconstruction is robust against motion and flow artifacts and can significantly reduce measurement time to 30 min at a resolution of 180 µm3 . Concisely, radial acquisition with combined compressed sensing and parallel imaging proved to be an excellent method for analyzing complex flow patterns in mice.

Histological Correlates of Diffusion-Weighted Magnetic Resonance Microscopy in a Mouse Model of Mesial Temporal Lobe Epilepsy.

Mesial temporal lobe epilepsy (MTLE) is the most common type of focal epilepsy. It is frequently associated with abnormal MRI findings, which are caused by underlying cellular, structural, and chemical changes at the micro-scale. In the current study, it is investigated to which extent these alterations correspond to imaging features detected by high resolution magnetic resonance imaging in the intrahippocampal kainate mouse model of MTLE. Fixed hippocampal and whole-brain sections of mouse brain tissue from nine animals under physiological and chronically epileptic conditions were examined using structural and diffusion-weighted MRI. Microstructural details were investigated based on a direct comparison with immunohistochemical analyses of the same specimen. Within the hippocampal formation, diffusion streamlines could be visualized corresponding to dendrites of CA1 pyramidal cells and granule cells, as well as mossy fibers and Schaffer collaterals. Statistically significant changes in diffusivities, fractional anisotropy, and diffusion orientations could be detected in tissue samples from chronically epileptic animals compared to healthy controls, corresponding to microstructural alterations (degeneration of pyramidal cells, dispersion of the granule cell layer, and sprouting of mossy fibers). The diffusion parameters were significantly correlated with histologically determined cell densities. These findings demonstrate that high-resolution diffusion-weighted MRI can resolve subtle microstructural changes in epileptic hippocampal tissue corresponding to histopathological features in MTLE.

Time-domain principal component reconstruction (tPCR): A more efficient and stable iterative reconstruction framework for non-Cartesian functional MRI.

PURPOSE: To improve the reconstruction efficiency (i.e., computational load) and stability of iterative reconstruction for non-Cartesian fMRI when using high undersampling rates and/or in the presence of strong off-resonance effects. THEORY AND METHODS: The magnetic resonance encephalography (MREG) sequence with 3D non-Cartesian trajectory and 0.1s repetition time (TR) was applied to acquire fMRI datasets. Different from a conventional time-point-by-time-point sequential reconstruction (SR), the proposed time-domain principal component reconstruction (tPCR) performs three steps: (1) decomposing the k-t-space fMRI datasets into time-domain principal component space using singular value decomposition, (2) reconstructing each principal component with redistributed computation power according to their weights, and (3) combining the reconstructed principal components back to image-t-space. The comparison of reconstruction accuracy was performed by simulation experiments and then verified in real fMRI data. RESULTS: The simulation experiments showed that the proposed tPCR was able to significantly reduce reconstruction errors, and subsequent functional activation errors, relative to SR at identical computational cost. Alternatively, at fixed reconstruction accuracy, computation time was greatly reduced. The improved performance was particularly obvious for L1-norm nonlinear reconstructions relative to L2-norm linear reconstructions and robust to different regularization strength, undersampling rates, and off-resonance effects intensity. By examining activation maps, tPCR was also found to give similar improvements in real fMRI experiments. CONCLUSION: The proposed proof-of-concept tPCR framework could improve (1) the reconstruction efficiency of iterative reconstruction, and (2) the reconstruction stability especially for nonlinear reconstructions. As a practical consideration, the improved reconstruction speed promotes the application of highly undersampled non-Cartesian fast fMRI.

Design of a shim coil array matched to the human brain anatomy.

PURPOSE: The purpose of this study is to introduce a novel design method of a shim coil array specifically optimized for whole brain shimming and to compare the performance of the resulting coils to conventional spherical harmonic shimming. METHODS: The proposed design approach is based on the stream function method and singular value decomposition. Eighty-four field maps from 12 volunteers measured in seven different head positions were used during the design process. The cross validation technique was applied to find an optimal number of coil elements in the array. Additional 42 field maps from 6 further volunteers were used for an independent validation. A bootstrapping technique was used to estimate the required population size to achieve a stable coil design. RESULTS: Shimming using 12 and 24 coil elements outperforms fourth- and fifth-order spherical harmonic shimming for all measured field maps, respectively. Coil elements show novel coil layouts compared to the conventional spherical harmonic coils and existing multi-coils. Both leave-one-out and independent validation demonstrate the generalization ability of the designed arrays. The bootstrapping analysis predicts that field maps from approximately 140 subjects need to be acquired to arrive at a stable design. CONCLUSIONS: The results demonstrate the validity of the proposed method to design a shim coil array matched to the human brain anatomy, which naturally satisfies the laws of electrodynamics. The design method may also be applied to develop new shim coil arrays matched to other human organs.

Direct modelling of gradient artifacts for EEG-fMRI denoising and motion tracking.

OBJECTIVE: Simultaneous electroencephalography and functional magnetic resonance imaging recording (EEG-fMRI) has been widely used in neuroscientific and clinical research. The artifacts in the recorded EEG resulting from rapidly switching magnetic field gradients are usually corrected by average-artifact subtraction (AAS) due to their repetitive nature. But the performance of AAS is often disrupted by altered artifact waveforms across epochs, notably due to head motion. APPROACH: Here, a method is proposed to make use of the known MR sequence gradient waveforms for a direct modelling of gradient artifacts. After accounting for filtering effects on the gradient artifacts, a continuous modulation of the gradient waveforms superimposed on the EEG signal is obtained. MAIN RESULTS: Although a moving AAS template can adjust to slow drifts in gradient artifact variation, it fails to adapt to abrupt motion, resulting in residual noise. We demonstrate how this modelling approach can reduce motion-affected gradient artifacts without distorting the underlying neuronal signals. Moreover, the method provides useful head motion information highly correlated with motion tracked by an optical camera. SIGNIFICANCE: Our work provides a novel way to improve gradient artifact removal in EEG-fMRI, and shows a potential to detect head motion without requiring additional hardware.